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1. Introduction 



In the theoretical limit of very large charm and bottom masses, Aqcd ^ ^ ^b, 
approximate symmetries of QCD allow one to express all six form factors relevant to 
semileptonic decays of B mesons into D and D* mesons in terms of the S-number form 
factor of the B meson, F{q^): 



where q = p — p' . S-number conservation implies F{0) = 1. Thus the six form factors for 
are known, in this limit, at one kinematic point. This is remarkable because 
it allows in principle for the extraction of the mixing angle Vcb from experiment without 
the need for a phenomenological model of the hadrons involved. In practice, however, the 
differential decay rate vanishes at the kinematic point where the form factors are predicted, 
so an extrapolation is necessary. The extracted value of the mixing angle depends on the 
nature of the extrapolation. Different models suggest different functional dependence of 
-F(g^), and we are stuck again with a model dependent extraction of the mixing angle. 

De Rafael and Taron^ have recently argued that the form factor F[q^) must satisfy 
some rather restrictive model-independent conditions in the form of bounds: 



where F±[q^) are known functions which depend on -F(O). If these bounds were correct, 
they would rule out some models which have been used to extract Vcb from experiment, 
making the determination of this mixing angle far more precise. 

Our discussion is organized as follows. In section 2 we review the derivation of Ref. |@]. 
In section 3 we note that their conclusion conflicts with the known physical behavior when 
applied to some other systems. Section 4 identifies the omission of Ref. 0: part of the 
states contributing to the analytic structure of F{q^) were neglected. We note why this 
omission is more severe for heavy-heavy quark states. We show how the bound may be 
modified to include additional contributions. The resulting bound is useful, however, only 
if one can estimate or compute the couplings of additional states. We explain how this 
works in models, and consider in particular the case of the 't Hooft model. In the last 
section we suggest some avenues of research that may yield somewhat model dependent 
but useful bounds. 



{B{p')\ln^h\B{p)) = F{q^){p + p') 



(1.1) 



F_{q^)<F{q^)<F+{q^) 



(1.2) 
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2. Derivation of the Bound 



Our derivation of ttie bounds follows closely that of de Rafael and Taron, except that 
we avoid strictly taking the limit nib — > oo and we work with arbitrary dimensions of 
spacetime, D. Consider the two point function of the heavy current = by^b 

iq^^q. - q^g^^.Mq^) = ^ j d^'xe^'^^ {TJ^{x)MO)). (2.1) 
In QCD, for L> < 6 it satisfies a once subtracted dispersion relation: 



1 , Imn(t) , , 



The absorptive part, Imn(q'^), corresponds to all the contributions to the right-hand side 
of Eq. ( |2.1| ) with a real (on shell) state between the two current insertions. This is a 
sum of positive definite terms, so one can obtain strict inequalities by neglecting some 
contributions to the sum. Concentrating on the term with intermediate states of BB 
pairs, one finds 

^-n*"') S (4.)(^-i./^r(£±i) (^{y^ - Ml), (2.3) 



where the function F is the analytic continuation of the form factor in Eq. (|1.1| ) for spacelike 
momentum transfer. The dispersion relation can be used to extend this inequality into the 
timelike region: 

^«?^) >- (4.)..-../^r(^) Lt t3/.(. : i^wi ■ 

For large 6-quark mass the two point function can be reliably computed from perturbative 
QCD on the timelike region, and in particular at = 0. With Nc colors, it is easy to 
verify that, to leading order in a^, 

xiQ') = (^;^r(3 - D/2) dx x^l - x)^ml + x{l - x)QT^'-\ (2.5) 

Setting (5 = 0, changing variables to y = t/4M^ and recalling that Mb /mi, = 1 + 
C^(^qcd/w{,), one has, to leading order in Aqcd/w^, 



OO 



dyk{y)\F\^ <1, (2.6) 
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where we have introduced 

kiy) = Cy-y'iy-l)^^-'y' (2.7) 

and 

2^/2+4iv,r(^)r(6^)' ^ ^ ^ 

Next we argue that any function F satisfying (|2.6|) and analytic except for a cut 
on the real axis for y > 1 is bounded as in ( |1.2| ), and we will explicitly construct the 
functions F± for any D < 6. To this end we map the complex y-plane onto the unit disk 
< 1 by the transformation 

V^ = ^^ (2.9) 

1 — z 

The cut y > 1 is mapped into the unit circle z = e*^. (The two edges of the cut are 
mapped into the upper and lower semicircles) The timelike region y < is mapped into 
the segment of the real axis < z < 1. In terms of this new variable the inequality ( ^.61 ) is 

^j\ew{e)\F\^<i, (2.10) 

where 

w{e) = k{y{e))^ 

de (2.11) 

= Ccos^ e/2 

As we will see shortly, the derivation of the inequality uses a function (p{z) analytic in 
\z\ < 1 such that |0(e*^)p = w{9). The construction of such a function is not unique 
— one may always multiply one such function by a power of z. We make the choice that 
reproduces the result of de Rafael and Taron: 

(Piz) = C'/^2-'/\l + ^)^/2(l - ;2)(5-^)/2. (2.12) 

With ref. [Q, let us define an inner product on the space of complex functions of a 
real variable 9, with < ^ < 27r, by 

if. 9)^^ r der{e)g{e). (2.13) 

Next let 

/i(^)=<^(e^^)F(e^^) 

UO) - (2.14) 

W) = 1 
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With this, we have 



1 

2^ 



2n 



ddw{d)\F\^ < 



1 



TT 



(2.15) 



Using Cauchy's theorem we can evaluate the other inner products. For example, with a 
contour C given by the unit circle \z\ = 1, 

27r 



1 

27ri 



dz 

c z- Zq 



{f2jl) 



(2.16) 







1 - zoe- 

Prom the positivity of the inner product we have that the 3x3 matrix (/i, fj) has positive 
determinant. This gives 

/ 9{zo) 9(0)' 
det I 9*izo) 1 I > 0. (2.17) 

g*{0) 1 1 



The inequality (|T3) 



follows, with 



F±{z) 



9(0) 



1 ± 



1 



n\g{0)\'^ 



The result can be used to obtain a bound on the charge radius, 



dF 

(.^) = 2(1) - 1)^(0) < 



D-1 

8M2 



^5 / 25 
^-2 + V^-^ 



(2.18) 



^2.19) 



De Rafael and Taron phrase their result as a bound on the Isgur-Wise function, by 
taking the infinite mass limit literally. We refrain from taking the limit nib = oo since 
for vv' > 1 this sends the form factor to zero. The Isgur-Wise function is not a physical 
observable, e.g., it depends on the renormalization point /i. In fact, the form factor F{q'^) 
is related to the Isgur-Wise function $,{vv',fx) through 

aL {vv') 

i{vv\^) (2.20) 



in the limit mh —>■ oo. In this limit only the leading log term needs be retained. In 
Eq. (|2.20|) the anomalous dimension is 



ttLivv') 



33- 2n/ 

where rif = 4 is the number of active fiavors and 

1 



vv'riyv') — 1] 



r{x) = 



Ini X -\- y — 1 



Notice that, for vv' > 1, the factor {asirrif,) /as{fx))°''^^'"'" ^ vanishes as mi, 



(2.21) 



(2.22) 



oo. 
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3. Trouble on the horizon 



De Rafael and Taron point out that some models of the form factor violate the bound 
( p. .21) ; they view this as a shortcoming of such models. 

One can turn this sort of observation around to cast doubt on the validity of the bound 
by trying to apply it first to simpler physical cases than the intricate heavy-light systems 
which are of primary interest. In particular, one may argue that the bound Eq. ( |2.19| ) 
cannot be correct because the charge radius of hydrogen atom-like wave-functions for the 
B meson must become arbitrarily large in the weakly bound case of the small coupling 
limit. 

Isgur[H has pointed out how to refine this objection further. The derivation above 
holds not just for the form factor of the B meson, but also for that of any particle that 
can be pair produced from the vacuum by the current J^. Thus, the bound must hold 
for the form factor of the Be meson, and even for that of the r]b- Now, these are well 
known to be well described by weakly bound non-relativistic wave- functions. Moreover, 
the charge radius is enormous compared to the Compton wavelength of the 6-quark, so that 
the violation of the bound is rather dramatic (as opposed to the rather marginal violation 
in the case of heavy-light systems). 

Clearly there is trouble, and a satisfactory explanation of what is wrong with the 
bound should account for the different behavior of form factors of heavy-light and heavy- 
heavy systems. 



4. The demise of the bound. 

Let us analyze in detail the assumptions and manipulations that lead to the form 
factor bound Eq. (|1.2|) . A number of physical considerations lead to the intermediate 
inequality ( |2.6D , which plays a crucial role. Is this inequality valid? As far as we can 
tell there is nothing wrong with the derivation that leads to it. In fact, the inequality 
is rather weak because one expects the absorptive part of the two point function, Imll, 
to be dominated by resonant intermediate states (T, T', etc). Is the perturbative QCD 
calculation of x reliable at = rather than in the more familiar limit Q oo? We 
believe it is. The onset of the physical cut is at spacelike momentum and scales with the 
square of the large mass m^,. Since this is much larger than the hadronic scale one expects 
QCD to be a very good approximation for timelike (and lightlike) momentum. 
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Consider therefore the second half of the derivation (foUowing Eq. ( 2.10| )) in which a 



positivity condition was derived for the function F{y). One of the key assumptions about 
the analytic structure of F — that F[q^) have no singularities below the S-threshold cut at 
> 4M^ — is in fact violated for good physical reason. This assumption is used to deduce 
(/2,/i) = g{zQ) in Eq. ( p.l6|) , and in computing (/3,/i) = g{Q)- Without this condition, 
the contours cannot be suitably deformed. We will argue that the form factor F[q^) has 
singularities for < 4M^, and that these modify the bounds in (1^) even in the large 
mh limit. 

The singularities we are referring to correspond to processes in which an 'onium' state 
is created by the current out of the vacuum with amplitude and then couples to the 
BB pair with amplitude QnBB- For simplicity, let us take these states to be stable so they 
contribute poles rather than cuts to the form factor. (The modification replacing poles by 
cuts and residues by integrals over discontinuities is straightforward and described below.) 
Since several of these states may lie below threshold, F^q^) has a term of the form 



When the q^ plane is mapped into the circle \y\ < 1, these poles land in the interior of the 
circle and hence in the interior of the contour of Eq. ( |2.16| ) . Now the naive scaling of these 
quantities for large rub is a factor of ^/mb for each heavy meson in an amplitude, so 

fn ~ ml^^ and g^sB ~ ^l^"^ , (4-2) 



and of course Mn ~ rrib. Therefore the contribution ([4.1| ) can be rewritten as 

yfngnBB ^43^ 

where = (2Mb)1/2/^, g^^^ = (2Mb)3/2^^^^, = 4M|M2 and, as above, y = 
q^ /AM"^. Clearly this contribution cannot be ignored even in the large mass limit. 

When these poles are taken into account Eq. ( |2.16| ) is no longer valid. Instead they 
generate an additional contribution: 



(/2,/i) = ^(^0) +g{zo), 
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(4.4) 



The explicit form of g for the case of pole singularities is 



= E 



n 



Ml{Zn- - Zn+) 



) 



1 



(4.5) 



where Zn± are the roots of z'^M^ - 2zM'^ + 4z + M!^ = 0, with \zn- \ < 1. 

A new bound is still obtained in this fashion, but without a priori knowledge of the 
quantities fn, QuBB ^^'^ it is hardly very useful. To obtain interesting constraints, one 
needs to first estimate the size of these contributions and second, to work with a system 
in which they can hopefully be shown to be small. 

Figure 1 shows the z-plane with the disc l^l < 1 corresponding to the complex = 
4M^y plane as per Eq. ( |2.9|) . The heavy line depicts the cut from some branch point below 
threshold out to infinity. It is clear from the figure that the angular integral in Eq. ( p.l6| ) 
and the contour integral cannot be made to agree, since the contour cannot be taken to 
cross the cut. An appropriate choice of contour differs from the angular integral by an 
integral over the discontinuity across the cut for \z\ < 1. Hence, again, Eq. ( |2.16| ) must be 
replaced as in Eq. (|4.4|) , (/2, /i) = g{zo) + g{zo), where now, in more generality. 



Here zr stands for the image under ( |2.9| ) of the location of the first branch point, yr = 
My/4M^ (discontinuities associated with light particles, e.g., tttv, are present but decouple 
in the large limit). 

We can now understand why the putative bound Eq. ( |2.19| ) was more severely violated 
for heavy-heavy mesons than for heavy-light mesons. As above we can phrase the question 
in terms of how large the slope of the form factor at = 0, the charge radius, may be. In 
ref. and ref. it was shown that the anomalously large radius of a heavy-heavy QCD 
bound state arises from the presence of many closely spaced singularities with significant 
residues. It is because the effect of such states was missed that de Rafael and Taron 
obtained an overly stringent bound on the charge radius of the heavy-heavy states. 

Furthermore, if one of these quarks is taken to be increasingly light, some of these sin- 
gularities migrate above threshold until, for a heavy-light system, one is left with only a few 
resonances below threshold (recall that the T*^'"^ is already above threshold). Therefore, 
the violation of the naive bound Eq. (|1 . 2| ) is less severe in these cases. Finally, we note that 
resonances continue to migrate above threshold for light-light systems, so the objections 
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Im 



(4.6) 



do not apply to the arguments of ref. (although the applicability of perturbative QCD 
in that case may cast a doubt on those results). 

Finally, we have examined all of these considerations in a relevant solvable toy model. 
We investigated the validity of these bounds in 1 + 1 dimensions in the large- A^c limit. This 
is an ideal testing ground for these issues since, on the one hand the quantities QnEB 
indeed the full singularity structure of the form factors can be computed exactly; and on 
the other hand much of the dynamics does in fact parallel 3+1 dimensional heavy quark 
physics in regimes where they can be compared. It is well known that this is a good 
approximation for Nc > 4. The bound in eq. (|1 . 2| ) , now with D = 2, holds for heavy- 
light systems provided Nc > 4, and is violated for Nc < 3. (This can be seen in Fig. 2 of 
ref. [§.) The bound is violated for heavy- heavy systems, even for rather large values of N^ 
recovering the observation above (but now without recourse to non-relativistic hadronic 
models) . 

In this model H] the absorptive part is saturated by resonances, so it has the form of 
an infinite sum over single poles on the positive axis, so that 

The coefficients and the masses can be calculated numerically. We have checked, 
numerically, that for > the perturbative expression, Eq. ( |2.5| ) with D = 2, is in 
excellent agreement with the sum in Eq. ( [4.7| ). 



5. Conclusions 

We have shown the form factor bounds on F[q^) for heavy-light mesons claimed by 
de Rafael and Taron do not hold because their derivation neglects contributions to the 
form factor below the S-threshold, for < 4M^. 

While the model-independent bound on form factors is lost, we can still derive a bound 
on form factors if we are willing to use a model to estimate the discontinuity across the 
cut in F[q^) for < 4M^. For example, a good approximation may be a sum of poles 
corresponding to the upsilon resonances below threshold, as in Eq. ( [4.3|) . If one can either 
measure or estimate the constants fn and g^BB that determine the residues, then a bound 
on the form factor is obtained. By including the effect of these very resonances on the 
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right hand side of ( |2.3|) an interestingly tight bound may foUow. This new 'model' of the 
form factor seems worth pursuing. 

It is also worth pointing out that nontrivial information about the form factors for 
the decay of heavy-light mesons into light-light mesons may possibly be obtained through 
similar methods. The argument is similar to that presented above, but now one considers 
the two-point function of a heavy-light current, = q'j^Q, and its hermitian conjugate. 
The disadvantage in this case is that one no longer has knowledge, a priori, of the value of 
the form factor at one value of q^. On the other hand, there is only one resonance below 
threshold, < q^^.^^^ = (Mb + M^)^. Thus, if one is willing to use, say, the constituent 
quark mo dell to calculate both F at Q'^ax = {Mb — Myr)^, and the two constants /s* and 
gB*B-Ki then one may obtain a bound on the form factor for all momentum in the physical 
region for the decay B neu. With any luck this bound will be tight enough that the 
form factor is determined to good accuracy. 
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^ Alternatively, some or all of these quantities may eventually be obtained from Monte Carlo 
simulations of lattice QCD, but a computation of the form factor over this wide momentum range 
may prove difficult on the lattice. 
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Figure Captions 

Fig. 1. 2;-plane with tlie disc \z\ < 1 corresponding to the complex = ^MqU plane 
as per Eq. (|2.9|) . The heavy Une depicts the cut from some branch point below 
threshold out to infinity. 
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